Enhanced Optical Dichroism of Graphene Nanoribbons 
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The optical conductivity of graphene nanoribbons is analytical and exactly derived. It is shown 
that the absence of translational invariance along the transverse direction allows considerable intra- 
band absorption in a narrow frequency window that varies with the ribbon width, and lies in the 
THz band for ribbons 10-100 nm wide. In this region the anisotropy in the optical conductivity can 
be as high as two orders of magnitude, which renders the medium dichroic, and allows near 100% 
polarizability with just a single layer of graphene. 

PACS numbers: 81.05.ue,81.05.ue,72.80.Vp,78.67.Wj 



Introduction 

Dichroism refers to the ability of some materials to 
absorb light differently, depending on the polarization 
state of the incoming wave, and leads to effects such as 
the rotation of the plane of polarization of light trans- 
mitted through them [1]. This characteristic is the ba- 
sis of several elementary optical elements like polarizers, 
wave retarders, etc., which are essential building blocks 
in photo-electronics and telecommunications. Dichroism 
is also widely relevant for substance characterization in 
fields ranging from spectroscopy, to chemistry, to life sci- 
ences. 

A grid of parallelly aligned metallic wires is a well 
known textbook example of a dichroic system, where 
unpolarized radiation becomes polarized perpendicularly 
to the wires, if the wavelength is much larger than the 
wire separation [2]. This example shows how geometri- 
cal anisotropy can be engineered to induce dichroism in 
otherwise isotropic media. 

Here we unveil the dichroic properties of graphene 
nanoribbons (GNR), and assess how effective grids of 
GNRs can be as polarizing elements. The motivation 
to explore GNRs in this context comes from a conver- 
gence of several critical properties. First, the optical ab- 
sorption spectrum of pristine graphene is roughly con- 
stant over an enormous band of frequencies [3, 4], from 
the THz to the near UV. This opens the unprecedented 
prospect of exploring its optical response to develop op- 
tical elements that can operate predictably and consis- 
tently in such broad frequency bands. Broadband polar- 
izers, for example, are a much needed element in photonic 
circuits for telecommunications, and graphene can play 
here an important role [5]. Second, the optical absorp- 
tion of graphene is easily switched on and off by varying 
the electronic density, which can be easily achieved by 
electrostatic gating [6] . Third, due to the record breaking 
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stiffness of the crystal lattice, one can suspend a graphene 
sheet and cut a grating of the thinnest nanowires (cur- 
rently of the order of 10 nm [7]), which opens new av- 
enues in ultra-narrow gratings, and upon which we base 
the system depicted in Fig. 1. Fourth, since graphene 
is metallic and a pure surface, the rich phenomenology 
associated with surface plasmons/polaritons is certainly 
unavoidable, further broadening the horizon of possibili- 
ties for optical applications. Finally, the atomic thickness 
of graphene results in a transparency of 97.7%. Hence, 
even if one is able to induce strong absorption along one 
direction, the overall transmissivity will still be large, 
which is important to maintain losses under control. 



Method 

The natural first step towards such possibilities con- 
sists in analyzing the intrinsic optical response of GNRs, 
to which we dedicate the remainder of this paper. We 
are interested in how the finite transverse dimension af- 
fects the optical absorption spectrum at low frequen- 
cies (IR and below), which is rather featureless in bulk 
graphene[3], but turns out to be richer in nanoribbons. 
The situation we envisage is depicted in Fig. 1, and con- 
sists in passing an electromagnetic wave across a grid 
of GNRs. For definiteness and technical simplicity we 
restrict our analysis to armchair (AC) nanoribbons, al- 
though our results do not depend on the specific chiral- 
ity, as will be clear later. Another important aspect to 
consider in GNRs has to do with how large edge dis- 
order is expected to be, and to which extent it might 
mask the phenomena under discussion. To address this, 
while at the same time keeping as much analytical con- 
trol over the results as possible, our calculations involve 
two steps. First, the frequency-dependent conductivity 
tensor crij{uj) of an AC GNR is derived exactly for free 
electrons governed by a nearest-neighbor tight-binding 
Hamiltonian (see below). The only energy scale is the 
hopping amplitude, t ~ 2.7 eV [8]. We then perform en- 
semble averages of such aij{uj), where the width is the 
fiuctuating parameter, and thus extract the overall re- 




FIG. 1: Illustration of the device geometry under consideration, consisting of a grid of parallel GNRs perpendicular to the 
incoming wave. The grid can be in vacuum, at the interface between two different dielectric media, or even inside a metallic 
waveguide with sectional area a x 6. A plane-polarized incoming wave has its polarization rotated by an angle 9 upon crossing 
the nanoribbon grating or, alternatively, an unpolarized wave emerges linearly polarized. 



sponse of the system accounting for "disorder". This 
procedure hinges on the assumption that the leading im- 
pact of disorder in the optical response is captured by 
the broadening of the quasi ID electronic bands, which 
is also achieved with an ensemble average of ribbons with 
fluctuating width. Moreover, such ensemble averaging is 
also close to the experimental situation, insofar as even 
state-of-the-art fabrication cannot control ribbon widths 
with atomic precision [9]. Thus, an array of ribbons cut 
out of a graphene sheet will always display a distribution 
of widths around a predefined target value {W) — Wq. 
Technically, the conductivity of such an array of GNRs is 
given by ((7,,(w)) = Y.w f{W)u'^{u:), where f{W) is the 
normal distribution for the ribbon width W , and a^ {ui) 
is the conductivity of a single ribbon of width W . 



Overall parametrizations are as follows. Ribbons are 
interchangeably characterized by their absolute width W , 
or by N , which counts the number of dimer rows along 
the transverse direction, and W = \/2>{N — l)a/2 ~ 
0.12 iVnm, where a ~ 1.42 A represents the C-C dis- 
tance. For the purposes of ensemble averaging, ribbon 
widths are uniformly distributed with a standard devi- 
ation that we take as constant, (iV^ — (iV)^)^/^ = 10. 
All the calculations discussed below have been done at 
T = 300 K. We use the terms intra- or mter-band in ref- 
erence to transitions occurring among mini-bands with 
the same or opposite sign of energy, respectively. The 
hopping amplitude sets the energy scale, and so we put 
t = 1 throughout. Conductivities are normalized to the 
universal value ctq — ne^ /2h of clean 2D graphene at low 
energies, and the incoming radiation has a wavelength 
much larger than W . 



Anisotropic Optical Absorption 

Lateral confinement, reduces the energy spectrum of 
GNRs to a set of mini-bands, each reflecting the disper- 
sion of an effective ID mode i {£ — 1,2,..., N), propa- 
gating longitudinally with momentum q: Eg^q^x = Xteg^q, 
where A = ±1, deflnes the valence and conduction mini- 
bands. 



Q,g 



= ^1 + 4 cos kg cos(g/2) -t- 4 cos^ kg 



(1) 



and ki = 7r^/(A^-|-l). Consequently, the density of states 
is dominated by Van Hove singularities (VHS) that de- 
velop at g = in the longitudinal momentum for each 
mini-band [10-12]. Such sharp spectral features trans- 
late in strong optical absorption for ideal GNRs, but are 
readily smoothed out by edge or bulk disorder and/or 
temperature in real systems. Our ensemble averaging 
has the same effect. 

In Fig. 2 we show the averages {(Jxx) and {(Jyy) for an 
ensemble with {N) = 150, and finite chemical potential 
/i = O.lt. As a consequence of time reversal symmetry, 
only the diagonal components of aij are non-zero. Trans- 
lation invariance along the longitudinal direction means 
that only inter-band transitions contribute to (Txx{'^) (see 
below). Consequently, {axx{'-^)) reproduces the bulk 2D 
behavior, as is clearly seen in the Figure. For the analy- 
sis of {(jyy{ijj)) is convenient to isolate the inter- (tr'"*™) 
and mtra-band (cr!,y*"*) contributions (the latter is al- 
lowed since along the transverse direction the electron 
scatters off the ribbon edges). Whereas (cr^'J,*™) feature- 
lessly follows {axx) (and hence the bulk 2D behavior), its 
mira-band counterpart displays a rather strong feature 
at low energies which, for this specific example, nearly 
reaches 14 times the universal value (Tq- 

Some aspects of Fig. 2 are worth underlying. Firstly, 
it is evident that, despite averaging to the same step- 
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FIG. 2: The three non-zero contributions for {au{ijj)), cal- 
culated for (N) = 150 (~ 18.5 nm), {N^ - {NfY^^ = 10, 
at r = 300A: and /i = O.lt (~ 0.3 eV). The inter-band 
contributions essentially follow the bulk 2D behavior, with 
a temperature-broadened step onset at fko = 2/i. In contrast, 
the mira-band contribution for the transverse conductivity is 
strongly peaked at low energies. The inset shows the same 
for a single ribbon of A'^ = 150. 



wise w-dependence, (o'y'^'°''(w)) is much smoother than 
(cr™*°''(a;)), even though the averages are over the same 
ensemble. This can be traced to the fact that, for each 
individual ribbon, only TV symmetric transitions {—E — > 
+E) contribute to {a'^^"' (oj)) , whereas (4"j,*°'(a;)) in- 
cludes 0{N'^) transitions among almost all pairs of mini- 
bands. Consequently, the latter has many more absorp- 
tion singularities, but much weaker, by conservation of 
spectral weight (this is explicitly shown in the inset of 
Fig. 2). The averaging is thus more efficient in wash- 
ing out the structure of VHSs in (o-J^'Jj^°''(w)). Secondly, 
the low-energy peak in (cr^'Jj^''^(w)) can be already identi- 
fied from a single ribbon (inset). Its origin is simple to 
understand with respect to Fig. 3. Since the bandstruc- 
ture consists of a set of discrete mini-bands, the chemi- 
cal potential will always be straddled by two of them at 
q — 0, such that £'f,g=o,A < M < -£'£+1,9=0, a- Given that 
transitions i —>■ £ + 1 are allowed in cr'"''''', one expects 
an absorption peak at fiw « |£'£,9,a — -E-^+i.g.Al- More- 
over, as per Eq. (10) the matrix element decays rapidly 
with the difference in band index, so that the transitions 
between the two bands closest to /i completely domi- 
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From Fig. 3 it is clear that there are al- 
ways two pairs of such bands, whose energy difference 
at g = is hSuji^2 ~ tta/B — ^^ ± 2fj./{N + 1). Since 
we are interested in situations where fi ^ t, the intra- 
band peaks are solely determined by the ribbon geom- 
etry: /iWniax ~ Tr\/3/{N + 1). This can be confirmed 
in the inset of Fig. 4 for ensembles with different (N), 
and introduces an element of predictability and tunabil- 
ity with respect to the frequency band where the optical 
absorption is highly enhanced. In other words, given the 
frequency of operation desired for a given application. 




FIG. 3: (a) Illustration of the transverse spectrum quantiza- 
tion and highlight of the two most significant inira-band tran- 
sitions 5a;i,2, which appear around the Dirac point (A*' — 21, 
Ep — OAt). (b) The projection of the mini-bands Ei^q^^ 
shown in (a), as a function of the ID momentum g, and for 
/ = 11, . . . , 18 (see text). 



one can select the appropriate ribbon width that yields 
the strongest optical anisotropy at the target frequency. 
To assess the polarizing efficiency of a single graphene 
ribbon we calculate the optical transmission amplitude, 
which is the ratio of the electric field amplitudes of the 
incoming and transmitted fields: ta{(jj) — Ea /Ea', 
{a = x,y). For radiation impinging normally upon an 
ensemble of GNRs separating medium 1 and medium 2 
(Fig. I), the transmission amplitude reads explicitly 



ta{u>) 



2Z(2) 



Z(1)+Z(2)[I + Z(1)(7„„(C^)]' 



{a^x,y), (2) 



where Z ~ -^//io/i/eoe is the impedance of each medium. 
Knowledge of ta{uj) allows the calculation of the degree 
of polarization, or polarizability V{uj), or the rotation of 
the plane of linear polarization {9 = 9f — 9i): 
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tan6'i, (3) 



The quantity V{uj) is particularly useful since we can 
immediately identify the degree of dichroism by how close 
|7^(a;)| is to unity (i.e. how close to an ideal polarizer are 
we). 

In Fig. 4 we plot the polarizability 'P(a;) for different 
ribbon widths and chemical potentials. It can be clearly 
seen that polarizabilities in excess of 50% can be achieved 
already with ribbons 45 nm wide. We underline that this 
is the polarizability of an atomically thin ensemble of rib- 
bons, which makes the magnitude of the effect even more 
striking! Even though the transparency of infinite 2D 
graphene is as large as 97.7%, the confinement-induced 
anisotropy can be so large as to suppress one of the field 
projections. The same figure also confirms that the op- 
timum polarizability is achieved at a width-dependent 
frequency Wmax which, as discussed above, has a simple 
form (inset of Fig. 4). However it is also clear that this 
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FIG. 4: The polarizability ■P(aj) in the low energy region, for 
different values of (TV) and fi ^ O.lt, T = 300K. The inset 
shows the position of the most prominent peak in crj^'y "^^ as a 
function of (A'') and /i. The n dependence is expectedly weak, 
while the peak position is seen to follow the analytical form 
described in the text. 



tunability is at the expense of the absolute amount of po- 
larizability (narrower ribbons — ^ larger ajmax ^ smaller 
7'(i^max))- Nevertheless, it has been experimentally con- 
firmed that the optical absorption of A^-layer graphene is 
trivially additive, from the bilayer to graphite [3, 13, 14]. 
This means that the effect reported here can be easily 
magnified by replacing single by few layer graphene, or 
by a succession of monolayer gratings. 

In addition, the form of Eq. (2) given in terms of the 
impedance of the media suggests that additional freedom 
can be achieved if the wave propagates inside a metallic 
waveguide. As is well known, electromagnetic propaga- 
tion in waveguides is restricted to normal TEM, TM or 
TE modes. Each of the latter two has a characteris- 
tic dispersion that is different from the free-space rela- 
tion CO = ck/n. For the purpose of analyzing transmis- 
sion and refiection amplitudes in a situation as depicted 
in Fig. 1, the effect of the waveguide can be absorbed 
in a renormalized and frequency-dependent impedance, 
Z{ijj). For example the mode TE^n has a character- 
istic impedance [15] Zm„(w) = Zu)/ y^oj'^ — uj^n^ where 
^mn = (c^7r^/A'e)[(TO/a)^ -I- (n/b)'^]. Hence each mode 
can only propagate if uj is beyond the mode cutoff fre- 
quency uJmm and this is frequently used to select/restrict 
the propagating modes by adapting the geometry of the 
waveguide. In our example we could take a square 
cross section (a = b), in which case the two degener- 
ate modes TEio and TEqi can be combined into an ar- 
bitrary incoming plane polarization [16]. In that case, 
if wio < cj < wii, only the modes TEio,oi propagate 

The 
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in the waveguide, and Zio{uj) — Zui/y/u)'^ — w^ 
cavity setup is interesting and useful for two reasons, 
which can be understood by inspection of Fig. 5: (i) on 
the one hand, by tuning the cavity dimensions so that 
"^lo ^ i^maxC-^) one can precisely cut the polarizability 
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FIG. 5: The effect of a metallic waveguide of square cross- 
section in the polarizability 'P{ijj) for two ensembles of ribbons 
((iV) = 150, 750). Each panel shows 'P{uj) for an incoming 
wave made of a combination [16] of the two lowest degenerate 
modes TEio,oi, with different cutoff frequency ojio. Each cjio 
is marked by a dot at the corresponding cj in the horizontal 



off below Wio, creating a well defined band of frequen- 
cies where the system displays high polarizability; (ii) 
on the other hand, since Ziq(uj) > Z (and, in particu- 
lar Ziq{uj > wio) 3> Z), the cavity highly magnifies the 
polarizability, even for a monolayer system. Taking as 
illustration the ribbon ensemble with (N) — 750 shown 
in Fig. 5, not only the polarizability V{uj) is made flatter 
and sharper, but its magnitude is considerably amplified 
in comparison with the freely-propagating wave {'P{lij) 
climbs beyond 80% in the entire frequency window). 



Derivation of the Conductivity 



The exact derivation of the conductivity tensor of 
an armchair graphene ribbon starts with the consid- 
eration of the nearest neighbor tight-binding Hamil- 
tonian describing the tt bands of graphene, and 
characterized by a hopping amplitude t ~ 2.7 eV 
[8]. The ribbon eigenstates have the analytical form 
[11, 17] |*,,,,a) = AA^„,„e-*9('"+"/2)sin(fc,n) x 
[\A,n,m) + Xe^'^'^'-i \B,n,m)), where kg, = t:1/{N+1) 
is the quantum number associated with transverse quan- 
tization (£ = 1, 2, . . . , TV), A/" = l/v^FTT, A = ±1 de- 
fines the valence (A =-1) or conduction (A =+1) bands, 
I A, 71, m) is the Wannier state at sub-lattice A of the unit 
cell at position R = nn + m m, N is the number of unit 
cells along the finite n direction, and q is the dimension- 
less momentum along m, whose value is within the range 
— TT < q < n. The phase difference between sub-lattices 
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arctan 



2 cos ki sin {q/2) 
1-1-2 cos kg cos {q/2) ' 



(4) 



This is sufScient to determine the optical conductivity 
from Kubo's formula [18]: 
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where 5' is the area of the ribbon, f{x) the Fermi dis- 
tribution function, and {^ i^q^\\ v^ \^ f.\q,\') is the matrix 
element of the /i component of the velocity (or dipolar) 
operator [19]. Since the energy scale is determined by i, 
let us introduce n — Hlj /t. 

Translation invariance along the longitudinal direction 
dictates that the matrix elements of the velocity Vx arc 
diagonal in q and £, leading to axx of the form 

3ff^ = AA, V 5^,,„AfJ(go, 4), (6) 

where 5fq,,^i,, = /{Ei^^qg^^) - f{Ee„,q,+), Ux = 
4/3-\/3(-/V — 1), and q^ is given by 



qa — 2 arccos 



{n/2y 



1 — 4 cos^ kf 



4 cos kf 



(7) 



The sum in Eq. (6) is restricted to those values of io such 
that qo e K. Finally, M^{qo,io) reads 



Mj(go,4) 



cos 9f, 



'■i^eo.qo-Qo/'^) cos ki„] 



sin(go/2) cosfcf 



(8) 



Only mfer-band absorption (i.e., transitions from the 
sub-bands with A = — ltoA = +l) contributes to axx- 

The analytical expression for ayy is slightly more cum- 
bersome than the previous one, due to the absence of 
translation invariance along that direction. As a re- 
sult, (i) the matrix elements of the operator Vy are non- 
diagonal in the sub-band index £, and (ii) there are both 
mfra-band (A = A') and mter-band (A ^ A') contribu- 
tions to the transverse conductivity. The calculation is, 
nevertheless, straightforward, yielding 
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^i^=^y2^2^ T^^ui2 StqXi. <(«0' ^1' ^2), (9) 

° f 1,^2 A, A' 

where Uy = 4/^/3(iV + 1){N^ - 1), ^^/^jj ,^ = 

nF{Ei,,q,,x) - npiEi^^q^x'), and Ve,,g, = 1 - (-1)^'+^^ 
This latter factor entails the selection rule for transitions 
among sub-bands ^i -1- ^2 = odd. The last factor is 
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sin kf, sin ke 



'yv-'— / sin2[(fc,^ +A;,J/2]sin2[(A:,, -fc,J/2] 
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Icosfc^^ef^,,;, + XX' cos ki^ei^^q J 



where C, 

go = 2 arccos 



?o/i/2 - 1 + XX' cos{ee^,q„ + 6i.^^qg - qo), and 
(02 - ai)Qb + VL^{bi + 62) ± Qc 



[hi - 62)2 



(11) 



with Q, = 2yfVL%ib2+n'^QbQa, Qb = h- 62, Qa = 

bia2 — 6201, fli = 1 + 4: cos^ ki- and hi = Acoski.. The 
sum in Eq. (9) is also restricted to those ^1,^2 such that 
qa € M, and to A < A' (photon absorption only). 

The expressions in Eqs. (6) and (9) are our central 
result, and from them follow all the averages and other 
physical quantities described and analyzed above. 



Discussion 

The optical absorption of a ribbon is seen here to 
be highly anisotropic on account of the new mtra-band 
channel made possible by the finite transverse direction, 
and the resulting electron scattering at the ribbon edges. 
Naturally, in the wide limit TV — >■ 00, the peak of ayy 
in Fig. 2 simultaneously narrows and migrates towards 
a; = 0, where it becomes the Drude singularity that we 
expect for an infinite and disorder-free system. Indeed, 
the easiest way to understand the sharp feature of ayy at 
low energies is to see it as a Drude peak that has been 
shifted to finite w by making the system finite along the 
transverse direction, thus allowing mira-band transitions 
of finite frequency. 

The issue of how to actually manufacture a grid of 
narrow GNRs with consistent and predictable width has 
been addressed earlier: it can be achieved by means of 
high precision patterning using a He-ion beam micro- 
scope in lithography mode [7], or more standard etch 
masks able to cut down to the 10 nm scale [20]. An al- 
ternative to cutting ribbons out of graphene sheets is the 
recently developed technique of unzipping carbon nan- 
otubes (CNTs) [21-23]. Nowadays is possible to pro- 
duce batches of CNTs with similar radius [24], and so 
this would allow the production of high quality ribbons 
without edge disorder. Another alternative, that com- 
pletely bypasses patterning, consists in inducing effective 
nanoribbons by engineering a non-homogeneous distribu- 
tion of strain in a graphene sheet [25] . 

As always in the context of GNRs, the role of disor- 
der needs to be addressed, and perhaps electron-electron 
interactions as well [18]. It is known that disorder can af- 
fect and even destroy many intrinsic features, such as the 
edge modes in zig-zag (ZZ) GNRs [10], the spontaneous 
spin polarization expected for ideal ZZ ribbons [12, 26], 
the width scaling of the gap [27, 28], or their conduc- 
tance [29]. In our case, disorder can modify the intrinsic 
optical anisotropy in difi'erent ways: (i) generalized dis- 
order will broaden the mini-band levels — which in turn 
broadens the intra-band absorption peaks — and will 
also induce a Drude feature in the conductivity around 
a; = 0; (ii) adsorbates and other impurities can intro- 
duce spurious absorption features; (iii) edge disorder can 
lead to localization of some electronic states [29]. With 
respect to (iii), much depends on the fabrication tech- 
nique, and the CNT unzipping method (or perhaps the 
strain-engineering route) would be preferred to mitigate 
edge disorder. Regarding (ii) post-patterning annealing 



techniques have been progressively improved, and proven 
quite efhcient in removing such sources of disorder [30]. 
Finally, the effects of generalized disorder (i) are expected 
to be less important for narrower ribbons, simply be- 
cause the anisotropy is induced by mtra-mini-band ab- 
sorption, and the separation of the mini-bands grows as 
oc 1/N . Consequently, unless disorder is strong enough 
to broaden the mini-bands by the same magnitude as the 
level spacing, well defined transitions are still expected, 
and the anisotropy should remain. Likewise, the appear- 
ance of a disorder-induced Drude feature at oj « is not 
expected to drastically affect the absorption peaks which 
occur at w = Wmax = finite. Moreover, as we pointed 
out earlier, by considering the response of an ensemble of 
GNRs with fluctuating widths, we are introducing con- 
siderable broadening effects already [compare, for exam- 
ple, the peak in {ayy{io)) for an ensemble in Fig. 2, with 
the flve times more intense peak of a single ribbon (in- 
set)]. For these reasons, we believe that the dichroism of 
GNRs remains considerably enhanced in the presence of 
realistic moderate disorder. 

It is worth highlighting the fact that, since the dichro- 
ism stems here from purely spectral considerations, the 
chirality of the GNRs is immaterial. In fact, all ribbons 
have the same scaling of the spectral features with TV, 
irrespective of their chirality, and so we can even expect 
the dichroism to remain when the ensemble comprises 
GNRs of arbitrary chirality. 

Finally, having in mind the depiction in Fig. 1 where 
we propose a grating of GNRs, we point out that this 
dichroism is intrinsic to each element of the grating. This 
is a departure from the conventional situation where the 
grating is made from a normal (isotropic) metal, and the 
polarizing effect arises from the geometry only, not from 
some intrinsic anisotropy of the metallic comb itself. In 
fact, it might have been noted that, whereas a conven- 
tional metallic grating polarizes perpendicularly to the 
slit direction, the dichroism of the individual GNRs fa- 
vors polarization along the ribbon direction. The ac- 
tual overall polarizing characteristics of a grating based 
on GNRs would have to be determined by the combi- 



nation of this intrinsic dichroism with the geometrical 
effect (just as in a conventional grating), and for which 
the surface plasmon-polariton physics plays an important 
role [31]. 



Conclusions 

Having derived the exact optical conductivity of 
GNRs, we studied the optical absorption response of 
ensembles of ribbons with fluctuating width. One ver- 
ifies that the optical absorption can be made highly 
anisotropic within a frequency band that is tunable via 
the ribbon average widths, and/or via the impedance 
characteristics of the embedding medium. Quantitative 
analysis reveals that an ensemble of monolayer GNRs can 
show near 100% polarizability, which is quite remarkable 
given the ultimate thickness of the polarizing element. 

The current analysis focuses on the intrinsic absorp- 
tion anisotropy of GNRs, where disorder effects are mim- 
icked by the fiuctuating ribbon widths. We are currently 
exploring routes to study the infiuence of real disorder, 
and combining the intrinsic absorption response of GNRs 
with the geometric effects expected to arise from the of 
the GNR grating. 

Given the recent developments in precision patterning 
and growth of narrow GNRs, and given the technological 
interest in optical elements operating in the IR and THz 
bands, we trust these results can motivate further theo- 
retical and experimental investigation of GNRs and other 
graphene-derived structures towards such applications. 
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